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Abstract. The Hipparcos Intermediate Astrometric Data have been used lately to estimate the inclination of the orbital plane 
of candidate extrasolar planets. Whereas most of these investigations derive almost face-on orbits, we show that the astrometric 
data are seldom precise enough to undertake such studies and that the 'face-on' result might be just a spurious effect of the 
method. 
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1. Introduction 

Even before the release of the Hipparcos and Tycho Catalogues 
( ESA 1997 ), their intere st regarding the extr asolar planets 
was already pointed out (Perry man et al. 1997). For the past 



two years, some investigators have re-processed the Hipparcos 
Intermediate Astrometric Data (IAD) of diff erent classes of ob- 
jects. Though initially limited to binaries ([Sfaderhjelm 1999 ; 
[Halbwachs et al. 20QC ; Pourbaix & Jorissen 2000| ), the planet 
hunters have also got interested in such re-processings (Mazeh 
etal. 1999; |Gatewood et al. 2000] [Zucker & Mazeh 2000| ). A re- 
cent statistical study based on the IAD concluded that a signif- 
icant fraction of the orbits of extrasolar planets are seen almost 
face-on (Han et al. 2001), thus pushing up their mass estimates. 
According to that study, a substantial number of these objects 
would no longer be planets. 

We have re-analyzed the IAD of 46 stars with planetary 
candidates and found that the orbital model, with the spectro- 
scopic parameters assumed, seldom improves the astrometric 
fit significantly. We show that in most cases the small incli- 



nations found by Han et al. (2001 ) are just an artifact of the 



fitting procedure. Such a caveat is indeed related to the size of 
the orbit with respect to the precision of the astrometric mea- 
surements. Similar results, although wrong, would therefore be 
derived wit h classes of instru ments of the same given precision 
(e.g. MAP (|Gatewood 1987b). 
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2. Astrometric model and the choice of the 
parameters 

Regardless of its physical nature, the motion of any member of 
a binary system with respect to the common center of mass is 



described by the well-known relations (Binnendijk 1960) 



x = A(cos E - e) + F \fl - e 2 sin E, 
y = B(cos E - e) + G \fl - e 2 sin E 



(1) 
(2) 



where E is the eccentric anomaly, solution of Kepler's equa- 
tion for given e, P, and T. Seven orbital parameters are thus 
required: A, B, F, G, e, P, and T. The first four, the Thiele- 
Innes constants, are linked to the Campbell elements (a a , i, u>i, 
and ft) through 



A = a a (coswi cosfi — sinwi sinf2cosi), 

B = a a (coswi sinfi + sinwi cos 51 cos j), 

F = a a (— sincJi cosfi — cost^i sinficosi), 

G = a a (— sina^i sinfi + cosc^i cosficosz). 



(3) 
(4) 
(5) 
(6) 



There is a priori no mathematical reason to prefer the Campbell 
elements over the Thiele-Innes ones. However, the spectro- 
scopic orbit yields, among other quantities, u>i and the prod- 
uct ai sin i (a/ is here reckoned in linear units whereas a a 
is its angular counterpart). Therefore, it is convenient to use 
Campbell's elements even if gl obal optimization technique s are 
then required to minimize x 2 ( Pourbaix & Jorissen 200C ) be- 
cause the model is highly non-linear in terms of the remaining 
elements. 

Most of the sub-stellar objects have been detected through 
radial velocity surveys and, therefore, the spectroscopic orbit 
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of the primary is known. Coupled to the parallax (m), two pa- 
rameters thus remain unknown, namely i and f2. Indeed, a a can 
be written as 



3.35729138 IO^^iPv 7 ! -e 2 vu/smi, 



where K\ is the amplitude of the radial velocity curve reck- 
oned in m s _1 , P is the period in yr, m in mas and e is the 
eccentricity. 

The IAD were made available to allow some further pro- 
cessing of the Hipparcos observations, especially those of 
binaries. Regardless of the model used to derived the final 
Hipparcos results, the IAD are residuals with respect to the 
single star model. When the binary nature of a star was dis- 
covered after the release of the Catalogues, the IAD can there- 
fore b e re-processed with an or bital model (Halbwachs et al. 
2000; Pourbaix & Jorissen 200C ). Although they are applied to 
1 -dimensional data, the equations of the orbital perturbation are 
essentially those of the visual absolute orbit. They thus supply 
with the inclination i. 

It is very tempting to re-process the IAD of all the stars 
around which planets were lately discovered. Indeed, the mo- 
tion of these stars should reveal tiny astrometric perturbations 
which can now be explained by the planet orbiting them. Fitting 
the IAD should supply with the inclination and therefore with 
the mass of the companion. For instance, that technique was 
applie d to HP 10697 (HIP 815 9) to obtain the mass of a brown 
dwarf (Zucker & Mazeh 2000). The brown dwarf nature of the 



companion was derived from the rather low value of the incli- 
nation, thus leading to a quite large mass for the secondary. 

3. Inclination from a tiny orbital wobble 

Why does fitting i and yield almost face-on orbits whereas, 
from a naive reasoning, one would expect preferentially 
edge-on situations for sp ectroscopically detected companions? 
Halbwachs et al. (2000 ) barely mention that feature. Indeed, 
in order to validate their approach, they process the IAD of 
known single stars with the orbital model (with ui±, e, P, and T 
fixed) and investigate the behavior of a a . Instead of a a — (the 
real value for single stars), they derive its Rayleigh distribution 
and conclude that a a is slightly overestimated (when the same 
method is applied to genuine binaries). 

They also noticed that the estimated a a is directly related to 
the residuals of the coordinates of the star, i.e. about 1 mas for 
Hipparcos. Fitting an orbit to the IAD when a a sinz <C 1 mas 
always leads to small inclinations: the smaller the product, the 
smaller the inclination. This is only due to the relative precision 
of the instrument with respect to the magnitude of the orbital 
wobble. The criterion depends on a a smi, not on K±. Indeed, 
long-period astrometric binaries also characterized by small 
K\ would be detected with a much longer Hipparcos-type mis- 
sion because, mutatis mutandis, a a sin i would increase with P. 
The orbital periods of known sub-stellar candidates are simply 
too small for a a sin i to be large enough. 

What about the 'precision' of such inclinations? Let us 
have a look at the derivative of the Thiele-Innes constant A 
with respect to i (the reasoning is the same for B, F, and G). 
These derivatives are indeed used to build the Fisher matrix 



whose inverse is the covariance matrix of the fitted parameters 
(Bevington & Robinson 1992). 

dA da a 



(7) di 



di 



(cos oj\ cos fl — sin ioi sin f2 cos i) 



a a sin uj\ sin f2 sin i. 



(8) 



If a a is an independent parameter in the astrometric solution, 
the first term disappears and the second vanishes with i. In that 
case, the smaller the inclination, the larger its uncertainty. On 
the other hand, when Eq. (Q) is adopted as a constraint, 



da a 

-— = -a a coti. 
oi 



(9) 



Hence, the smaller the inclination, the larger the derivatives, 
the smaller the standard deviation of the inclination, regardless 
of the precision of the astrometric data. 

4. Do these orbits improve the fit? 

Fitting the IAD o f HP 209458 (HIP 108859) using the spec- 
troscopic orbit by Mazeh et al. (200C) yields i « 0.02° even 



if we know that transits do occur (Charbonneau et al. 2000), 
implying i ~ 90°. We need a way to rule out such solutions. 

Table [j] lists extrasolar planet candidates with their a a sin i 
based on their Hipparcos parallax and their most recent spectro- 
scopic orbit. For all but seven, the value of a a sin i constrains 
the inclination to be < 10° (or > 170°) if one assumes that 
Hipparcos noticed the orbital motion. That latter assumption is 
where the weakness of the reasoning lies. An F-test can be per- 
formed to see whether the two additional parameters improve 
the fit of the IAP with respect to the single star model. In the 
absence of an astrometric wobble, the quantity 



F = 



N-7 X 2 s-X 2 c 



(10) 



follows the F-distribution with (2, N — 7) degrees of freedom 
( Bevington & Robinson 1992 ). N denotes the number of data 
points and x% an ^ Xc Slle tne va l ue °f *he \ 2 w i m th e 5- 
parameter (single star) model and orbital model respectively. 

The column labeled a in Table |l]gives the probability of ob- 
taining an P-value greater or equal to F if the null hypothesis 
HO: 'no orbital wobble present in the IAP' holds. That hypoth- 
esis is rejected for all but four stars at a 5% level. Even if two 
rejections are expected by chance in this sample, HP 38529 
(HIP 27253), HP 83443 (HIP 47202), p CrB (HIP 78459), and 
HP 195019 (HIP 100970) do deserve some further investiga- 
tions. However, a small a only means that an orbital model 
improves the fit of the IAP, with no guaranty that the adopted 
parameters indeed yield the best possible fit. So the obtained 
inclination might be unreliable even when a rj 0. 

The F-test rejects the astrometric orbital solutions for v 
And (HIP 7513) and for HP 10697 as well as the face-on so- 
lution we would obtain for HP 209458. It is worth pointing 
out that the value of a depends on the adopted parameters. 
Therefore, an alternative orbit might still substantially improve 
the fit even if the value of a listed in Table [l] is high. 

One should also mention that a a sin i large with respect to 
the precision of the instrument is not a sufficient condition for a 
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Table 1. List of planet candidates and their a a sini based on the spectroscopic orbit and the Hipparcos parallaxes, i is the 
inclination derived from the IAD when a a sin i is used as a constraint, a is the probability of obtaining an F-value greater or 
equal to F if there is no orbital wobble present in the IAD. Ref is the reference for the orbit: 1: Naef et al. (2000); 2: Queloz et al. 
(2000a); 3: [Butler et al. (1999|); 4: |Vogt et al. (2000|); 5: [Fischer et al. (2001|); 6: |Queloz et al. (2000b[); 7: |Marcy et al. (2000|); 8: 



Kurster et al. (2000|); 9: |Hatzes et al. (2000[ ); 10: [Butler et al. (2000[); 1 1: [Butler et al. (1997[); 12: |Naef et al. (2001[); 13: M ayor 
et al. (2000); 14: |Korzennik et al. (2000^ 15: [Butler & Marcy (1996[ ); 16: [Marcy et al. (1999| ); 1 7: |VIarcy & Butler (1996| ); 18: 
Udryetal. (2000a|); 19:[Noyes et al. (1999| ); 20: |Udry et al. (2000b[); 21: pochran et al. (1997[ ); 22: |Sivan et al. (2000[ ); 23: Mazeh 
et al. (2000); 24: pelfosse et al. (1998[); 25: |Mayor & Queloz (1995[) 
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reliable astrometric solution. For instance, Pourbaix & Jorissen 
(2000) fitted the IAD of 81 single-lined spectroscopic bina- 
ries and obtained reliable results for only 24 of them. For 21 
out of these 24, a a sini ranges from 0.6 to 13.8 mas whereas 
only 20 among the 57 others are characterized by a a a sin i 
much smaller than 1 mas. However, 27% of the stars rejected 
although they fulfill the 1 mas criterion have periods exceeding 
1 1 years. So no orbital solution could be derived for them due 
to the poor orbit coverage during the Hipparcos mission. 

In the present sample, HD 168443c (HIP 89844) and e Eri 
(HIP 16537) both have a a sin i > 1 mas but fail to get a signif- 
icant astrometric solution probably because of their relatively 
long periods (4.57 and 6.86 years, respectively). 



5. Conclusion 

Fitting (i, il) to the IAD when the spectroscopically con- 
strained a a sin i is much smaller than the astrometric precision 
always yields low values of sin i, irrespective of the true incli- 
nation. Although one cannot rule out the possibility of almost 
face-on orbits, very few of these orbits result in a significant 
improvement of the astrometric fit with respect to the single 
star model. 



Hipparcos was a very successful mission and its files cer- 
tainly still hold some hidden results. However, if one does not 
take care, one may be tempted to make these observations tell 
more than what they actually can. Moreover, the fact that one 
derives the same result with two instruments belonging to the 
same class of precision, for instance MAP and Hipparcos, does 
not always suffice to assess its reliability. Before a 100/ias-class 
instrument becomes available, the astrometric techniques will 
not derive the mass of most of the extrasolar sub-stellar com- 
panions known today. The four stars with a formally significant 
astrometric orbit (a < 5% in Table [j]) do on the other hand de- 
serve serious further consideration. 
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